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Abstract This paper is devoted to discuss the difference in the thermodynamic entropy 
budget per baryon in each type of stellar object found in Universe. We track and discuss 
the actual decrease of the stored baryonic thermodynamic entropy from the most prim¬ 
itive molecular cloud up to the final fate of matter in the black holes, passing through 
evolved states of matter as found in white dwarfs and neutron stars. We then discuss the 
case of actual stars of different masses throughout their evolution, clarifying the role of 
virial equilibrium condition for the decrease of the entropy and related issues. Finally, 
we discuss how gravity ultimately drives composition, hence structural changes along the 
stellar evolution all the way until the ultimate collapse to black holes, which may increase 
dramatically their entropy because of the gravitational contribution itself. 
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1 INTRODUCTION 

Entropy, as defined by physicists, is a mathematical function that encodes the thermodynamic macro¬ 
state of a physical system constructed from the statistical description of that system, thus a thermody¬ 
namic potential from which all the physical quantities may be computed. The classical works of Gibbs 
and Boltzmann (see a review by Muller 2007) clarified the meaning of entropy and suggested that it is 
also related to the degree of disorder of a system, although only in a very restricted sense, then mainly 
related to the heat capacity analogies with solids, liquids and gases’. Entropy drives the thermodynamic 
evolution of a system in time; and it is related to the amount of energy, including heat, that is available 
to do work according to Clausius and other contemporaries (Muller 2007). See also a discussion about 
the minimization of energy and the maximization of entropy and the relation to the available free energy 
in Muller (2008). 

It is often stated that entropy plays a key role in any process in Universe, and by means of its study 
that we can achieve a better understanding of the fate of the Universe and its contents. The Universe 
contains several differentiated structures from the largest to the smallest scales, and it is precisely among 
the latter that we are going to study the relations of entropy in various evolutionary states of the most 
fundamental astrophysical objects in Nature: stars, beginning with molecular clouds all the way down 
to the (ultimate) formation of black holes. The features of the thermodynamic entropy of a star along its 
main evolutionary phases will be addressed in Section 2. We present and discuss the results of entropy 
calculations in Section 3. We conclude in Section 4 by discussing how the whole stellar evolution is 
seen from the point of view of entropy considerations. 

* Metals have lower heat capacities and lower entropies than liquids. At the same time, metals have less ways to spread out 
some injected energy through their internal structures which makes them more orderly than liquids. 
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2 BASIC HYPOTHESES AND MODELS 


Stars are essentially self-gravitating systems held by some internal pressure against collapse. The stabil¬ 
ity of stars is due to the fact that they spend almost all their lives in stationary states, in which the virial 
equilibrium relation Epot -f ^Ekin = 0 among the various forms of potential Epot and kinetic Ekin 
energy terms hold. During most of their “active” lives stars generate energy through nuclear reactions, 
either in its simplest form (hydrogen to helium conversion) or advanced versions (helium to carbon and 
beyond). A set of time-scales describing how the structural, thermal and energetic adjustments are made 
can be defined, including 


- 


SttGp 


1/2 

is the free-fall time-scale, an upper limit to the maximum velocity of propaga¬ 


tion of any perturbation; 

- Tth ~ where Dth is the thermal diffusion coefficient (ratio of the conductivity to thermal ca¬ 
pacity), cfiaracterizing the time the star takes to establish a stationary distribution of the temperature 


when the latter is perturbed; 

- tkh ~ is the Kelvin-Helmholtz time-scale, related to the time it takes for the star to radiate 
of a significant fraction of the available energy; 

- Tnuc ( -jF ^ ) related to the burning of a given nuclear fuel with mass abundance X. 


Stars in steady state satisfy Tjj < Tth < Whenever a nuclear fuel is exhausted, 

the last inequality is violated, and the star seeks a new equilibrium state by contracting on a Kelvin- 
Helmholtz time-scale. Note that, because of the high thermal content and the first inequality, this is not 
actually a collapse. However, in each of the stages gravity gets stronger, and ultimately drives changes 
in the state of matter inside (i.e. degeneracy). Therefore, we can state that entropy should be produced 
but also radiated away (to the envelope and later away from the star), while hydrostatic equilibrium is 
maintained. In the formation of a star, or late in the final true collapse (supernova stage), this entropy 
generation and radiation is even more marked, since the contraction is much less “gentle” and out-of¬ 
equilibrium, therefore irreversible processes play a major role. 

We are thus lead to consider the Second Law of Thermodynamics in a familiar form for the system 
star + environment, namely the equation 1 



( 1 ) 


In principle, tracking all the sources of entropy inside the star (S) and the flux of the entropy currents 
Js throughout the boundary H we could calculate the increase or decrease of entropy for each stage of 
evolution of a given star of a given mass. Instead of that, one can just calculate initial and final states, 
thus “weighting” the relative importance of both terms on the right hand side of the above equation. Note 
that there could be ejection of mass (and entropy with it) in some explosive stages, although we shall 
not discuss the details of this complication in the remaining of the discussion, since we will compare 
just the entropy content in the final configurations of the objects under study (these final configurations 
will have the same number of baryons, for the reasons explained later). Besides, as we shall see later in 
Figure 1, a considerable amount of entropy is lost with the processes that lead to the explosive stages. 
In addition, the luminosity equation for a differential shell to evolve whenever there is a compression or 
expansion reads (in Lagrangian coordinates) 


^ 

dm dt 


( 2 ) 


with e the nuclear energy generation rate and S the entropy per mass unit and we have neglected 
the energy loss in the form of neutrinos. As we see, there could be luminosity generation even without 
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nuclear reactions, provided there is enough temporal variation of the enclosed entropy. Processes like 
ionization driven by compression, for example, contribute to the second term in some stages in between 
steady burning stages. The second term of Equation 2, also known as the thermal term, is more important 
for giant branch stars, not discussed here. 

A couple of comments are necessary regarding the above equation since it bears a very important 
meaning when one has to compute the structure of a star. If one turns off the term for nuclear energy 
generation, the virial theorem guarantees that the star will contract somewhat releasing part of its grav¬ 
itational energy to compensate the decrease of its internal energy. However, even if one does not “turn 
off” the source term, nuclear reactions slowly change the composition and the temperature gradients 
inside the star and then the structure, leading it to a new gravitational configuration as the star expands 
or contracts (depending on the internal energy balance). 

These gravitational adjustments imply that a gravitational work is done on the stellar matter which, 
in turn (due to the energy sources) drives an exchange of heat between adjacent shells of stellar mat¬ 
ter. In this way, the above equation is a direct consequence of the principle of conservation of energy: 
dQ/dt = dU/dt + PdV / dt. Thus, the change of entropy with time is a consequence of the very process 
of evolution of the star, i.e., the attempt of gravity to sustain the star in a state of (quasi) hydrostatic equi¬ 
librium, a heat-exchange process among matter shells (that is why some authors call this term, not very 
precisely, the “gravitational energy source”). Therefore, this term is related to tkh, and since the stars 
evolve on a much slower time-scale than tkh, except for the star on the Hertzprung gap, TdS/dt ~ 0 
and the condition of (local) thermal equilibrium is satisfied. Any complete and realistic model including 
transient adjustments must take into account the TdS/dt term. However, when computing the structure 
of a star we realize that we have to make some assumptions about the initial and boundary conditions 
and, in the case of a collapsing or expanding phase, the models may depend on these conditions, some¬ 
times quite strongly. Here we meet the real physical meaning of the equation just described: one cannot 
compute the structure of a star without knowing its previous history as stated in Clayton (1968). 

As a matter of fact, it should be remembered that the above eq. 2 does not determine the luminosity 
of a star as seen by a distant observer. That equation is prescription of what the luminosity should be 
in order to maintain the energy balance, e.g., the nuclear fusion and gravitational adjustments balance 
the energy losses, and is valid in any differential shell inside the star . The energy outflow from the 
star is ultimately determined by the radiation transport mechanisms such as diffusion, convection and 
conduction, all dependent on the shape and value of the temperature gradient, to be determined as a 
solution of the full set of structure equations. 

After all these considerations we shall, in order to be able of compare the different evolutionary 
phases of our model stars, first set a conserved quantity. The baryon number is precisely tailored for 
such a purpose, and we shall fix it to the value = 1.61 x 10®^ throughout this work, unless explicitly 
stated. This is not an arbitrary choice: the number corresponds to a mass (for small binding energy) 
of about ^ 1.35 Mq. It corresponds rather well to a solar-type example while hydrogen burning is 
considered, and even beyond, and is very close to the critical transition mass between the “evolutionary 
phases” from the high-mass tail of the mass distribution of white dwarfs to the low-to-average mass tail 
of the mass distribution of neutron stars as well. 

Our program to discuss the whole changes of the entropy will be the following: taking into account 
that a star will end up as one of the three kinds of compact objects (a white dwarf, a neutron star or a 
black hole),” we set these final configurations as having the the same number of baryons (our conserved 
quantity), N = 1.61 x 10®^, and we calculate the thermodynamic entropy of each compact object. After 
that we track back what would be the progenitor of each of our compact objects, e.g., what is the mass 
of the main sequence star that produced the compact objects with 1.61 x 10^^ baryons. We calculate 
the entropy of these stars in some chosen epoch of their lives during the Main Sequence. And then, we 
track back the entropy of the primordial clouds that produced these stars in the main sequence. 


^ Because the fate of a star depends on the mass it has at the moment it enters in the main sequence, the progenitors of the 
compact objects must have different initial masses. Roughly, a white dwarf has a progenitor with 1 — 7.5 Mq, a neutron star has 
progenitor with 8 — 25 Mq and a black hole has progenitor with > 25 Mq. 
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Our assumptions on the stellar state at each stage are: 

- the total energy is given by Etot = Eint + Ekin + Epou 

- the virial condition Epot = —2 x Ekin is satisfied; 

- the components have equilibrium particle distributions, for example, equipartition of energy holds 
for ideal gases, Ekin ^ kT. This is a justified as long as the thermal timescale remains very short, 
as is usually the case. 

In general, the entropy is a function of the (internal) energy, the volume and the number of particles 
of a system: S = S{Ein, V, N). Thus, we must properly choose the physical models that yield the 
energy and volume in each stage for a fixed number of baryons. The evolutionary stages we will discuss 
are: 

1 . White dwarf: White dwarfs are the endpoint of the evolution of ordinary stars with ^ 1 to ^ 7 Mq. 
From the point of view of stellar evolution, things depart considerably from the previous stages 
of the lives of ordinary stars. At this stage hydrogen can no longer burn in the star and the core 
contracts under its own gravity after achieving the Schonberg-Chandrasekhar condition (Schonberg 
& Chandrasekhar 1942) changing the hydrostatic equilibrium. This is the end of the Main Sequence 
stage (see below) and the star moves out of it entering into a completely new path in the H-R 
diagram. For a low- to intermediate mass star like the one we are dealing with, the contraction of 
the core proceeds until the point it is eventually halted by the degeneracy pressure of the electron 
gas, while the conservation of energy plus the virial relation together force the envelope to expand. 
After a series of structural changes, including the ignition of helium in degenerate conditions in the 
core (helium flash), the star will eject the outer envelope in a series of thermal pulses of increasing 
amplitude. Mass loss at this stage is very large and can not be ignored, therefore to hold the baryon 
number fixed as before we are not considering the actual evolutionary path but rather an ideal model 
situation for the sake of clarity. This is the end point of the evolution of this star because the remnant 
cannot generate energy. It will cool down releasing all the thermal energy it had stored. 

It is widely known that the actual composition of a “typical” white dwarf is mainly carbon-oxygen 
(C-O) and they have typical masses of about 0.6 Mq. However, in order to conserve the baryon 
number of 1.6 x 10®^ we had to assume the creation of a very heavy white dwarf with ~ 1.35Mq 
well to the tail of the mass distributions of these stars. The composition is then different: it is 
probably oxygen-magnesium-neon (O-Mg-Ne) white dwarf. To produce such a heavy white dwarf, 
the progenitor star was chosen to have mass of 7 Mq, roughly 8.3 x 10^^ baryons, which lose 
6.7 X 10^^ baryons during the mass loss phase, which is a reasonable assumption. 

Thus, in order to calculate the structure of this star we assumed the ultra-relativistic regime with 
a polytrope of index n ^ 3, since the mass of this object is quite near the Chandrasekhar limiting 
mass (Chandrasekhar 1931). The equation of state is then given by P = The assumed 

central density is 1 x IQ^^g/cw? from which Rwd — O.OO2i?0 (^ ISOOfcm) is obtained. 

We calculated the entropy of the WD in two specific moments: the hot initial phase, after the thermal 
pulses, where the core temperature is Thwd ^ 5 x 10® A" and for a very late and evolved phase, 
when the core temperature is Tcwd ^ 1 x 10® A". 

2. Neutron star: Neutron stars are are formed by the collapse of a massive star 8 — 25Mq resulting 
in a compact object of > 1.2M0 and R ^ lOfcm®. Neutron stars are supposed to have all the same 
composition since the burning stages of the evolution reach the limit of the iron, from where no 
more exothermic processes are possible. So, with iron as the starting point for further evolution, the 
core contracts to a completely new phase of (dense) matter. 

Neutrons stars also have a mass distribution, possibly two peaked (Valentim et al. 2011), and the 
different masses are possibly due to the masses of the progenitor stars and the mass loss processes 
during earlier stages of evolution during the post main sequence phase. Thus, in order to produce a 
neutron star with 1.6 x 10®^ baryons, near to the lighter part of the mass distribution, we assumed 


^ Another possible way to form a neutron star is via the accretion-induced collapse of a white dwarf (van den Heuvel 2011). 
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a progenitor with 11 Mq that lose something about 11.5 x 10®^ baryons in the ultimate supernova 
explosion. 

We calculated the entropy of our neutron star in three distinct moments; the hot proto-neutron star 
phase whose temperature was assumed to be T ^ 5 x and Rpns ^ 55.75km (this is 

about five times the radius of the forthcoming neutron star, due entirely to the “hot phase”); a later 
“stationary” hot phase with temperature of T = 1 x 10® A" and Rns — 11.15/cm, which settles 
a few hours after the formation at most; and a final cold “stationary” phase with temperature of 
T = 1 X 10^ A" and Rns — 1115/cm representing the cooling of the isolated neutron star after 
^ lO^yr approximately. 

3. Black hole: We end our calculations with the ultimate state of the collapsed matter, the black hole. 

The actual formation of this extreme compact object is marked by the death of a very massive star 
(> 25Mq). In a similar way as with the neutron star formation, we assumed a progenitor with 
25Mq that, after its normal evolution, ejects 28.1 x 10®^ baryons ending with the formation of 
a 1.6 X 10®^ baryons black hole. After the formation of the event horizon, the final object emits 
thermal radiation at a Hawking temperature Tbh = 2 kbTrG(Nm ) ~ ^ lO^^AT. The celebrated 

proportionality between the entropy and the area (Bekenstein 1973, 1974) now applies, since all 
forms of matter have disappeared beyond the horizon. 

4. Main sequence star: We follow the entropy evolution, dominated by the ideal gas component, 
using Townsend’s MadStar online tool (http;//www.astro.wisc.edu/~townsend/) to create ordinary 
main sequence star models, burning hydrogen to helium and solar metallicity. We created four main 
sequence stars: one with 1.35 Mq that will evolve in some 4 billions of years into an old ordinary 
star (it will evolve further to a C-O white dwarf that we will NOT study here); one with 7 Mq that 
will evolve in some hundreds of millions of years to a O-Mg-Ne white dwarf; one with 11 Mq that 
will evolve in some dozens of millions of years to a neutron star; and one with 25 Mq that will 
evolve in some hundreds of thousands of years to a black hole. The MadStar is a online tool based 
upon an approach by Bill Paxton on the famous Eggleton code. Although some limitations exist, 
none of them have strong influence to our purposes. 

5. Molecular cloud: Modern determinations of stellar forming conditions (Caproni et al. 2000) have 

shown the occurrence of substantial clumping within molecular clouds. We consider the formation 
of stars inside these clumps which have typical temperatures ~ 20Ar, typical masses ~ 0.2Mq and 
typical densities ~ that merge together to form a single star. The radius of each small 

cloud is then Rsmc — 8.23 x 10^®cm ~ 0.03pc. Thus, each small cloud amounts to 1.19 x 
10®® baryons. The clumps at the immediate stage of star formation are opaque to radiation, but 
their entropy is largely dominated by the ideal gas component (we neglect magnetic fields in this 
discussion). The merging of a few of these small clumps will produce a star with the assumed 
baryon content, in a complex process driving the star towards the Zero-Age Main Sequence (ZAMS) 
immediately following the ignition of hydrogen and the establishment of the hydrostatic equilibrium 
condition. 

Recalling, we want to study, at the same time, the thermodynamic entropy of a given 1.6 x lO®’^ 
baryons in different states of matter, i.e., in different degrees of compactification and the changes 
of the entropy content along the the evolution of the objects that originated these final compacted 
baryons in the first place. As mentioned above, we have three compact stars representing three 
different exotic states of matter. We need then four original molecular clouds: three for the three 
compact stars and another one to account for 1.6 x 10®^ baryons enclosed in a ordinary star in 
“normal” state. 

Clumping of small molecular clouds add up to form four molecular clouds, with 1.6 x 10®^ baryons, 
8.3 X 10®^ baryons, 13.1 x 10®^ baryons and 29.7 x 10®^ baryons, respectively, that will form the 
four main sequence stars that, in turn, will form our final four objects whose entropy will be studied. 

It is important to check for each stage of evolution the state of degeneracy, since degenerate gases 
follow a different entropy expression than ideal gases. In Table 1 we give the Fermi temperature in each 
stage, remembering that the degeneracy occurs if Tobject << Tp. 
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Table 1 Fermi temperatures [K] of each evolutionary stage given by our hypotheses: molec¬ 
ular cloud (MC), ordinary Main Sequence (MS) star, white dwarf (WD), proto-neutron star 
(PNS), and neutron star (NS). It is important to notice that in the table below we calculated an 
average Fermi temperature employing the average density of each object. The Fermi temper- 

ature is given by Tp = 2 ^ (’ where rj = y, generally a function of the radial 

coordinate r. 



MC 

StarO 

Starl 

Star3 

Star4 

WD 

PNS 

NS 

1 + 

4 * 

~ 10 '' 

~ 5 X 10 "“ 

~ 10 '’ 

- 150 

~ 10 '’ 

- 120 

~ lO'" 

~ 70 

~ 5 X 10 '" 

- 30 

~ 10 "" 

~ 5 X 10 ’’ 

~ 10 "-" 

~ 10 "“ 

~ lO"'^ 

~ 10 "^ 


Employing the values in Table 1 we can calculate the entropies in all the stages. 


3 RESULTS AND DISCUSSION 

In each subsection below we discuss how we calculate the components of the entropy and its total value 
in each evolutionary stage. Then, in Table 2 we show which term of the entropy is dominant. This is 
of course impossible for a black hole, which does not have any ordinary component left and needs a 
separate consideration. 

A comparison with other known entropy sources in the Universe is interesting, as discussed by 
Frampton et al. (2009). It is important to remark that stars giving rise to neutron stars/black holes rep¬ 
resent ^ 1% of the 10^^ stars present in the visible universe. In any case, the entropy content of all the 
stellar populations is tiny compared to other known components (i.e. CMB photons) and thus irrelevant 
for the whole budget. 

3.1 Molecular cloud 

Giant molecular clouds are the main cradle of stars. Fragmentation of a giant cloud and further clumping 
and collapse of smaller units form the main blocks, as discussed above. These are composed mainly by 
neutral molecular hydrogen, being described roughly by the ideal gas law. Here we assume that the small 
clouds described in Section 2 clump together to amount the determined number of baryons as described 
in the previous section and that each is in equilibrium just before the clumping and collapse to form the 
stars in the zero-age main sequence (ZAMS). Because the temperature is Tmc = 20Ar ^ T'pe- jions 
there is no degeneracy. Then, the entropy can be calculated by the expression for an ideal gas: 


SmC = ^ SsMCbaryons = ^ NsMC,bh f ) ' 2 


5 ln( ^inSMC \ 

V NsMC,b ' 


■ const 


( 3 ) 


where const = ^ -I- |, k\, is the Boltzmann constant, Vsmc is the volume of each small 

cloud and niu is the atomic mass unit. The important assumption here is that the clouds are composed 
by molecular hydrogen only. 

The total entropy in each of our four cases is the sum of the entropy of a certain number of small 
clouds that add up to a specific number of baryons. So, for our first object, 6.73 small clouds add up 
1.6 X 10®^; for the second, 34.98 small clouds add up 8.3 x 10®^; for the third, 54.97 small clouds add 
up 13.1 X 10^^; and for the fourth, 124.93 small clouds add up 29.7 x 10®^. Recall that each small cloud 
has 1.19 X 10^® baryons, as described in the previous section. 

From our assumptions: Etot = 0, Epot = —2Ekin and Ekin = we finally find 

Q 

Smci = 4.77 X 10‘^‘^erg/ K or S ^ -—— 

kbN 


= 21.47; 
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and 


SmC2 

SmC3 


24.8 X lO'^'^erg/K or 
39.0 X lO^^erg/K or 


S ^ 


S 

S 

kbN 


21.65; 


21.57; 


SmC 4 = 88.6 X IQ'^^erg/K or S ^ = 21.62. 

kbN 

The calculation of the entropy of a molecular cloud in an earlier stage is tricky and requires a careful 
consideration (not attempted here). In the transparent stages of the cloud, the radiation is not effectively 
coupled to matter and it is not clear whether it should be included. Nevertheless, this stage happens well 
before any actual condensation stage and is not important for our considerations. 


3.2 Main Sequence stars 


The second evolutionary stage encompass the Main Sequence (MS) stars, the region in the Hertzprung- 
Russell (HR) diagram where the stars stay most of their lives. For stars in the range of masses worked 
out here, 1.35 Mq to 25 Mq, the period of residence in the MS is ^ while the energy 

generation is mainly due to the so-called p — p chain for star with masses up to 2 Mq and due to the 
CNO cycle for masses above 2 Mq. In our case, this is roughly 4 Gyr. As stated, typical temperatures 
of our four models are > 10^K while the Fermi temperature is ^ 10^K for electrons and ^ 10^ iT for 
ions of hydrogen Therefore, there is essentially no degeneracy along the Main Sequence, except 

maybe for a small degree of degeneracy in the inner core. 

The four molecular clouds collapsed to form four models of Main Sequence stars with ZAMS 
masses 1.35 Mq, 7 Mq, 11 Mq, and 25 Mq. The MadStar evolutionary code calculates the entropy of 
the structure in a way that the entropy of each star at Main Sequence is: 


S'i.35Mq 


8.26 X 10'^*erg/K or 


S ^ 


S 


= 3718; 


'S'tMq = 53.5 X lO'^'^erg/K or 


S->- 


S 


= 4671; 


5'iiMq 


78.3 X m^'^erg/K or 


S ^ 


S 

kbN 


= 4331; 


S 25 Me = 212 X lO^^erg/K or S ^ = 5173. 

The general behaviour of entropy with ageing inside the MS is to get smaller and smaller, as illus¬ 
trated with the example of the star with 1.35 Mq: 

= 8.26 X lO^Vp/i^ or 5 ^ = 3718; 

kbN 


St=0.8 Gyr = 8.13 X lO'^'^erg/K or 


S ^ 


S 

kbN 


3659; 


5t=2 Gyr = 7.70 X 10^'^erg/K or 


S'-> 


5 

kbN 


3466; 


and 
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Gyr = 5.10 X 10^'^erg/K or S ^ = 2295. 

A more general consideration of the behaviour of the entropy inside the Main Sequence for several 
masses and in the Post-Main Sequence stages will be given in section 3.7. 

3.3 White dwarf 

The condition of matter inside a white dwarf is fairly different in the two stages we considered. In the 
hot phase, the electrons are degenerate, but the ions constitute basically a Boltzmann gas. Therefore we 
have three components for the total entropy. While we can still use the same terms in the expression 3 

for the ions and use the expression Srad = for rh® radiation, we need a new expression for 

the entropy of the degenerate matter; 

SwDe - 7 - 2 -’ 

I 

where Xe = and p/e- is the Fermi momentum of the electron sea. The factor 1/2 comes from 
the supposition that Np = Nn = We, that is reasonable for a white dwarf. From our calculation 5'e- ~ 
IQ'^^erg/K and Sions ^ lO'^^erp/iT. 

Thus, for the hot initial state we obtain 

SwDhot ^ 1.27 X IQ'^'^erg/K or S ^= 5.72. 

kbN 

However, as the star cools down the electrons remain degenerate, but the ions suffer a phase tran¬ 
sition to form a Coulomb lattice (Mestel & Ruderman 1967). The Debye temperature marking the 
crossover of these regimes is for our model 9d — 1-8 x lO^AT. Because we chose a final “station¬ 
ary” state with a temperature ^ 10^A" << 9jj, vte are well inside the regime where the heat capacity 
goes with ~ corresponding to phonon lattice excitations. The entropy is given then by 

/ \ 3 

- lO^^erg/K (5) 


= 1.86 X 10“®. 

kbN 

We see that the entropy budget changes along the cooling age of a white dwarf. At first, in the 
hot WD state, most of the entropy is stored in the ions. In the cold WD state, the degenerate electrons 
eventually hold the largest fraction of entropy, being the excess of entropy carried away by the photons. 

3.4 Neutron star 

Again, the state of matter inside the star differs radically in these phases. In the proto-neutron star 
phase the assumed temperature is 5 x 10^^ AT resulting in electrons which are still degenerate (the Fermi 
temperature in this configuration is Tpe- ~ lO^^W), but the neutrons can still be considered a non¬ 
degenerate gas with some degeneracy correction (the Fermi temperature in this state is Tpions ^ 5 x 


Sw Dcold—ions — 


IGNiT^kb / T 
15 I ^ 


while for the degenerate electrons S^- ~ IQ^^erg/K. 
For the cold final state we then obtain 


SwDcoid — 4.14 X lO'^^erp/AT or S 
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Then, the entropy of a proto-neutron" star is given by the entropy of degenerate electrons plus the 
entropy of a neutron Boltzmann gas; 

SpNS - 2.96 X or S ^ 7 ^- = 13.32. 

kbN 

After the proto-neutron star phase, the neutron star settles and quickly cools down via neutrino emis¬ 
sion. Our model for cold neutron stars results from solving the Tolman-Oppenheimer-Volkoff equations 
(Weber 1999) complemented by the SLy4 equation of state (Douchin & Haensel 2001), a particularly 
suitable choice for a compact star composed by very neutron rich matter with interactions. We selected 
the star with N = 1.6 x 10®^ baryons which correspond to a (gravitational) mass M ~ 1.23M0 and 
radius i? = 11.15fcm, showing the effects of stronger gravity in these objects through a lower total mass 
(larger binding). 

Assuming a temperature of T^s 1 x 10® AT in this later stage, the neutrons become degenerate and 
support the star against further collapse. Then the stellar entropy is given by the entropy of degenerate 
neutrons (we have neglected the small fraction of electrons/protons enforced by beta equilibrium): 

n^{xl + iy/^NkJ^) 

Sns = - 2 -’ 

where x„ = and pfn is the Fermi momentum of the neutron sea. 

The entropy is then 


SNShot ^ 1.21 X lO^^erg/K or S ^ = 5.45 x 10"^ 

kbN 

As the neutron star cools down, eventually reaching a temperature of Tjys 1 x 10^K in the core, 
its entropy decreases further to 

SnscoM- 1-21 X 10^^erg/K or ^ ^ ^ = 5.45 x 10“^ 

kbN 

At this stage, the entropy decrease is mainly due to photon emission. 

3.5 Black hole 

A black hole is a region in space-time in which an event horizon has been formed and enclosed all 
the matter of the progenitor, and ultimately it does not matter which kind of particles contributed to it. 
Because there is no access to the interior content of a black hole, a thermodynamic description of the 
collapse cannot be based on the entropy of the contents since these are lost from the observable universe. 

The black hole entropy depends solely on the observable properties of the black hole; mass, electric 
charge and angular momentum. Because of the area theorem (Bekenstein 1973, 1974), these three pa¬ 
rameters appear in a combination defining the area. The expression for the entropy of a Schwarzschild 
(non-rotating, uncharged) black hole is given by: 


Sbh = 


kbA ‘fkb'nG{NmuY 


AGh he 

where we used the fact that A = AirR'jj = 167r(G'M/c®)® with Rh = 2GM/c^. 

The numeric value of the entropy of 1.6 x 10®^ baryons enclosed by the event horizon is then 


(7) 


Sbh — 2.63 x ffr^erg/K or S 


kbN 


= 1.2 X 10 


,20 


4 


Here we used a radius Rpns = ^R-ns^ where Rns — 11-15/cm; which mimics the hot stage before cooling 
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That is many orders of magnitude larger than the preceding states of matter. This can be considered a 
measure of the degree of irreversibility of the ultimate collapse to form black holes out of known matter 
and energy, that is, the entropy currents are quickly stopped from flowing out the black hole because 
of the horizon formation and the highly irreversible collapse generates huge amounts of entropy. This 
qualitative picture, however, cannot be taken too seriously, since it is not clear yet how exactly the 
entropy is located in the area (making S a non-extensive quantity), a subject of much discussion and 
calculations for string theorists and loop quantum gravity researchers (Ghosh & Perez 2011; Jacobson 
et al. 2005) beyond the scope of our work. 

3.6 Thermodynamic entropy summary 

In Table 2 we show the dominant entropy contribution in each compact object with 1.6 x 10®^ baryons 
reflecting, at the same time, how entropy end up stored in the final configurations at different states 
of matter and the evolutionary path as the stars go through their lives. In Figure 1 we show the ther¬ 
modynamic entropy per baryon in units of kf, as a function of the central density at each stage of the 
evolution (notice that for a black hole we assumed an effective central density of 10^®^/cm® for plotting 
purposes). 

It is most important to realize that in stellar evolution gravity drives ultimately the changes in the en¬ 
ergetic processes in the interior of stars and the final states of matter in the compact objects accordingly 
to the mass of each progenitor star. 

The important comparison is then between the progenitor main sequence stars of 7 Mq, 11 Mq 
and 25 Mq, and the “initial” final configurations hot white dwarf, hot neutron star (passing through the 
proto-neutron star phase), and black hole. We clearly see a trend in lowering the entropy from earlier 
stages until the final configurations. 

From the point of view of the baryon content in relation to the state of matter inside the compact 
objects, we see that although the main sequence star of 7 Mq has lower entropy than the 11 Mq star, the 
entropy jump to a state of lower entropy is higher from the star of 11 Mq to the neutron star than from 
the star of 7 Mq to the white dwarf. Thus, for the same number of baryons, the more compact object has 
the lower entropy and this difference is due to the structural changes of matter, from ordinary matter in 
white dwarfs to very neutron-rich matter in neutron stars. In this respect the entropy of the proto-neutron 
star is very similar to the entropy of the hot white dwarf, since in broad terms, the proto-neutron stars is 
a kind of iron white dwarf, with matter in a state that resembles the state of white dwarf. 

Gravity is the ultimate force driving the entropy changes and stellar evolution. Then one can wonder 
why the most compact object, the black hole, has a entropy that is so many orders of magnitude higher 
than the other compact objects with the same number of baryons, apparently contradicting the very 
conclusion we stated in the previous paragraph. The answer lays in the gravity field itself and its putative 
entropy content. We shall see below that the latter could nicely explain the big difference of entropy 
content between low and high curvature stars. 


Table 2 Entropy components [erg/K] of each compact object with 1.6 x 10®^ baryons given 
by our hypotheses: white dwarf (WD), proto-neutron star (PNS), neutron star (NS), and black 
hole (BH). 



Radiation 

Ideal baryons 

Ideal electrons 

Degenerate e 

Crystal 

Degenerate n 

Area 

HWD 

- 10®'' 

TTup 

- 

TTup 

- 

- 

- 

CWD 

- 10®® 

_ 

_ 

~ 10®® 

- 10®® 

_ 

_ 

PNS 

~10"® 

~10^® 

- 

~10"® 

- 

- 

- 

HNS 

~ 10®® 

- 

- 

- 

- 

~ 10®® 

- 

CNS 

- 10®® 

_ 

_ 

_ 

_ 

~10®" 

_ 

BH 

~ 10-1" 

- 

- 

- 

- 

- 

~10®i 
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Fig. 1 Thermodynamic entropy per baryon in units of kb versus the central density of the 
objects in each stage of evolution. The (red) diamonds are the four molecular clouds; the 
(green) squares are the four main sequence stars (note that there are five squares, but one 
of them represents an evolved version of the same 1.35 Mq star); and the (pink) circles 
represents the compact stars (white dwarf, neutron star and black holes plus the proto-neutron 
star). The (black) arrows show the changes (or “evolution”) of the entropy as stellar evolution 
proceeds. The two (blue) arrows going from the cold white dwarf to the proto-neutron star 
and from the cold neutron star to the black hole are the special case of induced collapse. 
Notice that for a black hole we assumed an effective central density of IQ^^g/crn? for plotting 
purposes 


To finish the thermodynamic entropy summary, we show in Figure 2 regions in the plane T vi p 
which entropy regime is dominant in the typical density and temperature range corresponding to the 
models studied here. We also show the well-known dominant pressure regimes for comparison. 

We see that the entropy regimes are not directly connected to the pressure regimes, although they 
overlap in many situations. This can be seen in the specific case of a white dwarf case: while the electron 
degeneracy pressure is dominant through the whole range of densities of white dwarfs, the entropy of 
the degenerate electrons is dominant only in the regime of low temperatures and high densities (the 
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Fig. 2 The region above the hrst left line (R) is radiation pressure-dominated; the region 
between the hrst left line and the second line (G) is gas pressure-dominated; the pressure of 
non-relativistic degenerate electrons (NRDe“) is dominant in the region between the second 
line and the hrst vertical line; the region between the two vertical lines is dominated by the 
pressure of the relativistic degenerate electrons (RDe^); the region to the right of the second 
vertical line is dominated by the pressure of degenerate neutrons (Dn). The other two lines, 
the ones crossing the hrst vertical line, dehne a change in the regime of the specihc heats of 
the ions inside white dwarfs due to crystallization. Regarding entropy, in the darkest regions 
the entropy is dominated by the entropy of the Boltzmann gas; in the second darker region 
the entropy is dominated by the entropy of a hot lattice of ions; in the third darker region the 
entropy is dominated by the entropy of cold crystallized ions lattice; in the lightest grey region 
of the WD box the entropy is dominated by the entropy of degenerate electrons; and hnally, 
in the region of the NS box the entropy is dominated by the entropy of degenerate neutrons. 


lighter grey region of the WD box of Figure 2). This change is directly related to the evolution of the 
white dwarfs, more specihcally to the cooling and to the fact that there is no energy generation inside 
the core of these objects. A similar behaviour may occur for the degenerate neutron composition case, 
and it is worth further study. 

3.7 Actual Main Sequence stars and beyond 

The considerations made above about the entropy of an idealized path for a hxed baryon number were 
intended to show how this quantity evolves as successive stellar equilibrium states arise. However, it is 
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clear that actual stars do not complete the whole evolution, and therefore it is important to discuss how 
the entropy behaves for the whole range of stellar masses in actual stars. 

It is currently agreed that nuclear reactions start for the central conditions of stars with mass above 
O.OSMq (for solar composition). Besides, the stellar structure is stable up to a high mass limit of at least 
90Mq and possibly of ^ 130Mq. Within these two extremes, the stellar structure varies and the entropy 
characterizing the star is subject to a quite different behaviour. Generally speaking, it is well-known that 
a stellar core or shell may transport energy in the radiative mode (diffusion of photons driven by a 
small temperature gradient, the so-called local thermodynamic equilibrium condition, LTE) or, when 
the stability criterion formulated in terms of the size of the temperature gradient ^ is violated, in the 
convective mode, in which large-scale motions of the fluid itself are involved (Kippenhahn & Weigert 
1994). The latter regime is difficult to model since a variety of scales and complex phenomena (i.e. 
turbulence) play a crucial role. Therefore, simple approaches ignoring the multi-scale dimension in 
favour of a representative “bubble” are employed, the most popular being the mixing length theory 
which parametrizes a number of physical effects with a single dimensionless quantity a, expected to be 
0(1). Even the plain substitution of the true temperature gradient by the adiabatic temperature gradient 
^ I is many times employed for calculations based on the fact that the actual value of the gradient 
cannot be large on physical grounds, since convection is very efficient and there is little room for a large 
value build-up of temperature differences (see the textbooks in Carroll & Ostlie (2006); Kippenhahn & 
Weigert (1994)). 

The behaviour of entropy in both cases (diffusive vs. convective) is very different, and an additional 
condition, that of virial equilibrium, is determinant for its value, at least as long as the balance is dom¬ 
inated by the ideal gas component and not by degenerate electrons or radiation, as we shall see in a 
moment. 

The virial relation states that, in equilibrium, the gravitational and thermal energies satisfy | = 

2^NkT = 3j^^kT, where N is the total number of particles, /i is the mean molecular weight of the 
gas and mn is the hydrogen mass. If we assume a constant value of the density, a simple manipulation 
of the relation yields T ^ Gfimn jpf2/3pi/3^ jjj jjjg expression of the ideal gas shows that 

the entropy decreases with increasing temperature, in complete agreement with the Second Law when 
subject to the virial equilibrium condition. It is said sometimes that the star “digs an entropy hole” 
along its evolutionary path, but this behaviour ceases as long as some other component (i.e. degenerate 
electrons) dominates the entropy balance. 

The observation that large-scale motion of the fluid is the dominating transport mechanism beyond 
a certain value of the temperature gradient (the convective regime) leads to a variety of situations along 
the stellar evolution for a given mass, and also to different configurations for stars with different masses. 
But another important feature of the stellar interior will be relevant for an overview of the entropy. It 
is related to the condition in which nuclear energy is released by fusion reactions. Simple calculations 
show that around ~ 2Mq the main reactions from CNO catalytic cycle overcome the so-called p-p 
channel, and therefore the heavier stars burn hydrogen in a much more dramatic way, since the latter is 
very dependent on the temperature (in contrast with the former which is very mildly dependent on T). 

Thus, stars below this “great divide” threshold develop a steep gradient only at the outer layers, 
whereas above it the steep gradient is present in the core. Therefore, the structure of stars in the first 
regime (termed Lower Main Sequence) is radiative on the inside, and convective outside, all the way 
to the photosphere. Those in the Upper Main Sequence, in turn, become convective inside and keep a 
radiative envelope. 

If we consider stars of lower and higher masses away from the threshold, there are two important 
boundaries to be noticed in the mass parameter: stars below ~ O.SM© are totally convective as the 
result of the systematic inwards advance of the convective envelope from ^ 2Mq solar mass towards the 
lowest values. In addition, the convective core grows with mass at the other end, i.e., Mconvective/Mtotai 
grows from ^ 50% for a 20 Mq star up to 80% for a 100 Mq star. Those stars in between present 
both convective and radiative regions, as discussed above. 
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Given these features and considering the homogenization provided by convection, the entropy pro¬ 
file of any star residing in the Main Sequence can be depicted as in Figure 3. The first extreme (below 
~ 0.3Mq) corresponds to a constant value of the entropy, as showed in the upper left panel of Figure 3. 
More massive stars, of the solar type, and up to around ^ 2Mq are represented by the upper right panel 
of the same Figure, featuring a radiative core and a convective envelope. Still more massive ones, up to 
^ 20Mq or so, correspond to the profile shown in the lower left panel; and those stars above this value 
are again represented by a constant value of S (upper left panel). 


Convective 

Radiative+convective 

Convective+radiative 

SeMereS shells with 

active nuclear reactions 


0.25 0.5 0.75 1 1.25 0.25 0.5 0.75 1 1.25 

R R 


Fig. 3 Entropy regimes (arbitrary units). For plotting purposes we used the same scale for all 
plots. 


From the point of view of the entropy, similar behaviour occurs when stars leave the Main Sequence 
evolving towards the right of the HR diagram. Low-mass stars develop an increasingly massive helium 
core, mostly degenerate for solar-mass stars and below. The envelope expands and becomes convective 
inwards (the reverse of the Hayashi track in the stellar formation process...), and the entropy distribution 
resembles the one in Figure 3 (lower right panel). Note that degenerate cores, by their own nature, 
do not diminish their entropies with increasing temperature. In contrast, upper main sequence stars 
develop normal, non-degenerate cores and depending on their exact value, end their lives as white dwarfs 
(those up to around 8Mq) or ignite further nuclear reactions. The full sequence of available nuclear 
cycles (including true fusion reactions and photo-disintegration rearrangements, yielding net energy) is 
achieved for masses beyond ^ IOMq. For the range 8 — IOMq it is expected that a degenerate O — 
Mg — Ne core is formed, but the conditions would not be enough to go beyond the carbon cycle. Since 
when each of the combustion cycles are no longer possible, the inert core contracts and the reactions 
formerly at the center migrate to a shell around it, the so-called “onion structure” develops for the more 
massive stars, in particular, in those with M > IOMq which would complete the ignition of all available 
nuclear reactions. The entropy, in turn, adopts a constant value inside convective cores, but given that 
the virial condition must be maintained, its value decreases for each new cycle. Meanwhile, the radiative 
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shells adopt the entropy distribution growing from the inner to the outer edge of the shell in each case. 
Thus, the entropy distribution for any given stage of massive stars is qualitatively similar to that shown 
in the lower right panel of Fig. 3. This process of decreasing the entropy per baryon at the center is 
more pronounced for the lighter range of massive stars, since the relation T ^ Gumn j\j2/3pil3 

to ^ oc and the ideal gas entropy then increases with increasing mass, being much higher for the 
range 15 — 20Mq and beyond. Since the end of these stars as core collapse supernovae depend on the 
binding energy of the core, directly related to the entropy per baryon, the low entropy per baryon in the 
~ 10Mq range was considered as a favourable condition for the explosions to succeed, although many 
other factors seem to be at play in this process (Murphy et al. 2013). 

At last, a remark should be made regarding actual stars and systems of stars. Binary evolution 
channels are very important for the formation and evolution of compact objects, but introduce new 
complications to the entropy considerations which lie beyond the scope of the present approach. For 
comprehensive considerations of a possible common envelope phase we refer to Ge et al. (2010a,b). 

4 CONCLUSIONS 

We have seen that the gravitational compactification of stellar matter makes the thermodynamic entropy 
of matter to decrease from the main sequence stars to their correspondent compact stars in an ordered, 
monotonic sequence. This decrease of the thermodynamic entropy of matter is not at all in contradic¬ 
tion with the Second Law since the gravitational contraction releases high-entropy radiation (and even 
neutrinos) in a way that more than compensates the decrease in the matter entropy of the object. Thus, 
it is the ability to radiate away entropy what keeps the entropy of stars decreasing. 

Starting with the condensations in a molecular cloud, as the collapses goes on, the core of the object 
is heated by the contraction up to the point the fusion reactions of protons into helium begins. At this 
moment, the processes in the core balance gravity and the contraction halts, but the thermodynamic 
entropy of matter remains lower than that it was in its initial state. Later on, as the nuclear fuel of the 
star is exhausted, a new contraction begins, being halted only when the heat capacity becomes positive 
again when a phase transition of matter occurs. After a few billion years, a white dwarf is formed with 
a still lower thermodynamic entropy. 

In other words, the entropy decrease of matter with the gravitational contraction is at the expense 
of an increase in the entropy of the immediate environment due to the release of very high-entropy 
radiation/neutrinos during the life of stars in the Main Sequence and beyond. There is a delicate balance 
of two processes. The fusion reactions per se can be thought to lower the entropy because the number of 
particles diminishes. However, this process is highly exothermic, and increases the local entropy due to 
local heating. Which effect dominates depends on the local temperature: at low temperature (relative to 
nuclear scales), the entropy gain from the exothermic reactions is favoured. When the temperature of the 
core reaches some threshold, the reactions become entropically unfavoured and stop. This coincides with 
the jump to the next stage of the evolution of the objects. In the white dwarf and neutron star stages, the 
entropy budget suffers a stronger influence from cooling processes, since the structure of these objects 
will remain unchanged for time-scales which are infinite in practice. However, the injection of entropy 
in the environment is never too extreme, for example, a (type II) massive star supernova will produce 
around 3 x lO^^erp jK of entropy mainly in neutrinos. A hundred times this figure will be injected due to 
the dissipation of the massive ejecta, ultimately leading to dust heating of the interstellar medium. These 
numbers are still lower than the total entropy introduced by the heating of dust due to the absorption of 
ordinary radiation of stars (Bousso et al. 2007; Frampton et al. 2009). 

Gravity plays a major role in all these contractions and entropy changes, driving stellar evolution 
towards its end. Gravity is responsible, again, for the formation of a black hole. As seen above, collapse 
under gravity causes entropy to increase enormously; however here the black hole itself “cools down” 
to a very low Hawking temperature (for a solar-mass or so black hole), much smaller than the external 
temperature while its entropy increases astonishingly. Because black holes have negative heat capacity, 
they absorb radiation faster than they can emit by the Hawking mechanism. We face a very special case: 
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the environment is cooled by the collapse (differently than in normal stars) while the entropy of the 
matter in the black hole is very high (Wallace 2010). 

As mentioned earlier, the gravitational field itself may carry and store entropy. This entropy has been 
related to the curvature inside the stars, but it is negligible from Newtonian objects like the molecular 
clouds, main sequence stars and white dwarfs. For neutron star, constructed under the framework of 
General Relativity, the entropy of gravity is expected to be larger, but still does not alter significantly 
the entropy budget of this very compact star. However, things change considerably when a massive 
main sequence star collapses to form a black hole. Black holes are “pure gravity”, a singular point in 
the metric that in which the curvature term goes to infinity, forming an object from which no matter 
can ever escape. That is why this collapse actually increases the entropy stored in this object. There is 
considerable activity in the community to construct and characterize a gravitational entropy. A definite 
recent example is the proposal of Clifton et al. (2013) in which a “gravitational entropy” has been 
suggested based on the Bel-Robinson tensor which makes use of the Weyl part of the curvature tensor 
Cabcd (see Sussman & Larena (2014) for applications). Thus, such a proposal attaches an increasing 
entropy to the gravitation which is related to the increasing curvature value. As a particular example, 
the gravitational entropy reduces to the Bekenstein-Hawking value for the case of a Schwartzchild black 
hole. Therefore, if this proposal stands, the origin of the black hole huge entropy could be thought as the 
limiting case of action of gravity in stars, not only by forcing the change of state of the matter (discussed 
above), but also leading to the highest available found in nature. 

The evolution of black hole in a cosmological environment has been considered before. At some 
point the background becomes as cool as the black hole itself, driven by the expansion of the Universe, 
but further evolution can not preserve the thermal equilibrium (Custodio & Horvath 2003). Rather, the 
black hole begins to evaporate and eventually disappears (Hawking 1975) or leaves a tiny quantum 
residual, meaning that most of its entropy returns eventually to the environment. 

We have seen in this work that entropy in its multifarious forms plays an important role in Stellar 
Evolution theory. We discussed how the gravitational contraction/collapse, although irreversible in na¬ 
ture lowers the entropy of matter relatively to its initial state. However, we still have a long journey 
towards complete understanding of the role of entropy in the fate of stars. We conjecture, however, that 
there is a synthesis to be made from the study of entropy in stars from a totally general point of view. 
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